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The transverse momentum dependent helicity distributions of valence quarks are calculated in 
the light-cone diquark model by adopting two different approaches. We use the model results 
to analyze the P^-dependent double spin asymmetries for 7r + , tt~ and 7r° productions in semi- 
inclusive deep inelastic scattering, and find that the asymmetries agree with the CLAS data in one 
of the approaches. By taking the Fourier transform of the transverse momentum dependent helicity 
distributions, we obtain the helicity distributions of valence quarks in the transverse coordinate 
space, and then apply them further to predict the Bessel-weighted double spin asymmetries of 
7r + , 7r~ and 7T° productions in semi- inclusive deep inelastic scattering at CLAS, COMPASS and 
HERMES for the first time. The shape of the Bessel-weighted double spin asymmetry thereby 
provides a direct probe on the transverse structure of longitudinally polarized quarks. 
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I. INTRODUCTION 

The quark helicity distribution, which measures the 
net helicity of quarks in a longitudinally polarized nu- 
cleon, 

9i( x ) = f q ym (x) - f q i/m (x) (1) 

is an essential observable encoding the spin structure of 
the nucleon [l|, 0]. During the last two decades, the x- 
dependence of the helicity distribution has been studied 
extensively, including the parametrizations [3|— rTj l from 
experimental data and model calculations I8 H10| . How- 
ever, the transverse momentum dependence [ill [l2j of 
the helicity distribution still remains poorly known. In 
recent years, it is found that the quark transverse motion 
plays a crucial role p^| - [T9| through spin-orbit correlation 
in certain spin phenomena, such as the single spin asym- 
metries measured in semi-inclusive deep inelastic scatter- 
ing (SIDIS) and elsewhere (see Ref. [26| and refer- 
ences therein). On the other hand, the inclusion of quark 
intrinsic transverse momentum can provide much richer 
information about the nucleon structure [271 . |28| that can 
not be obtained from the collinear picture. Thus it is 
necessary to revisit the quark transverse motion effect 
to the quark helicity distribution, especially its trans- 
verse momentum dependence that can help to build a 
3-dimensional partonic structure of the longitudinal po- 
larized nucleon. 

As a step toward a better understanding of this is- 
sue, we study in this paper the transverse momentum 
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dependent (TMD) helicity distribution of the valence 
quarks, denoted by gf L (x,p^), using the light-cone di- 
quark model, which has been applied to calculate the x- 
dependence of the helicity distribution Q , the transver- 
sity distribution I29l . |3fJ and a number of other TMD 
distributions [3ll - [34| . We will calculate the TMD helicity 
distribution in two different approaches corresponding to 
two different options of the unpolarized TMD distribu- 
tion. The TMD helicity distribution enters the descrip- 
tion of SIDIS process when the transverse momentum 
of the final hadron Phi is measured. Using the model 
results from the light-cone diquark model, we calculate 
the double spin asymmetry (DSA) of pion production in 
SIDIS process as function of Phi ■ The CLAS Collabora- 
tion has measured [2~ij such asymmetry using a polarized 
electron beam to collide a longitudinally polarized proton 
target recently, and such measurement provides the first 
experimental glimpse of the TMD helicity distribution. 

Although it is possible to access the TMD helicity dis- 
tribution from SIDIS data, usually one has to model the 
quark transverse momentum dependence in phenomeno- 
logical analysis because calculating the cross section dif- 
ferentiated by Phi. involves complicated convolution of 
transverse momenta (35j . Conventionally, the convo- 
lution can be resolved through the P/^-weighted pro- 
cedure 

nana 

in a model-independent way, and the 
weighted cross section involves the transverse moments 
of the TMD distributions. However, much of the infor- 
mation about the explicit transverse momentum depen- 
dence is lost in this way. An alternative approach has 
been developed u sing the Bessel functions as weighting 
functions in Ref. [37|, where the authors show that the 
convolution becomes a product of TMD distributions and 
TMD fragmentation functions (FFs) in the Fourier space 
(also referred to as the transverse coordinate space, de- 
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noted by by) conjugate to pr- The TMD distributions 
in the b^-space thus can be measured directly by the 
Bessel-weighted asymmetries whose expressions in lead- 
ing twist have been given in Ref. [13. This motivate us 
to examine the Bessel-weighted DSA in SIDIS by model 
calculation since it might be measured easier than other 
Bessel-weighted asymmetries. To this end we calculate 
the Fourier-transformed TMD helicity distribution in the 
light-cone diquark model and make the prediction for 
Bessel-weighted DSA of pion production in SIDIS for the 
first time. The results are presented at the kinematics of 
CLAS, HERMES [H and COMPASS [H, where there 
are ongoing SIDIS programs or there are data that could 
be analyzed. 

The paper is organized as follows. In Section II, we 
briefly introduce the calculation of the TMD helicity dis- 
tributions of valence quarks in the light-cone diquark 
model and discuss their properties. We then calculate the 
DSAs of 7r + , 7r~ and ir° productions in SIDIS. In Section 
III, we calculate the valence helicity distributions in the 
Fourier space and present the prediction of the Bessel- 
weighted DSAs of three pions in SIDIS at CLAS, HER- 
MES and COMPASS. Finally we provide summaries in 
Section IV. 



II. TMD HELICITY DISTRIBUTIONS IN 
Phx-DEPENDENT DOUBLE SPIN 
ASYMMETRIES 

The TMD helicity distribution giL.(x, p\) is defined 
through the matrix elements of the bilocal operator 



9xl{x,P t ) = - 



2 J (2tt) 3 



(2) 



where £(0, £) is the gauge link ensuring the gauge invari- 
ance of the definition. In this calculation the gauge link is 
set to be trivial since the helicity distribution is T-even. 

One can insert a complete set of intermediate states 
{n} into Eq.© 

91l ( X ,pI) ~yj t^S^ ((1 



(27r) 3 2£„ 



x)P^ 



P H 



x 5 2 (p T + P nT ) {(PS L \V+^ i+) (0)\n) 2 
-(PS L \V^ {+) (0)\n) 2 } , (3) 

where ip(+) — ^7~7 + V' is the "good" component of the 
quark field, and V± = \{\ ±75) are the helicity projec- 
tion operators. Similarly the TMD unpolarized distribu- 
tion has the form 



fi{x,p 2 T ) =—= 



d 3 P n 



-5 ((1 - x)P+ P+) 



V2-W (27r 3 )2£„ 

v n \ / 

x 5 2 {p T + P nT ) {(P\V + ^ (+) (0)\n) 2 

+ (P\V-A+m\ n ) 2 } ■ (4) 



The TMD helicity distribution has been calculated by 
different models (IoWIsT ] and by lattice QC D Here, 
we adopt the light-cone diquark model [3, |45| to calcu- 
late the valence TMD helicity distributions. In the SU(6) 
quark-diquark model, the proton state with a spin com- 
ponent S z = +h in the instant form can be written as 



|Pt) 



1 



(ttd) V - v^ud) V - V2(uu) d^ 



-2(uu) 1 d J -] + cos9(p s (ud) s u t , 



(5) 



where the two quarks inside the brackets form the scalar 
(S) or the vector (V) diquarks, f/l denotes the posi- 
tive/negative helicity of quarks, 9 is the mixing angle 
that breaks the SU(6) symmetry when 9 ^ 7r/4, and 
ifD{x,Px) is the wave function in the momentum space 
for the diquark D = S or V. Thus in this model the di- 
quarks are treated as the intermediate states appearing 
in Eqs. (jSJ) and (gj). 

Notice that the proton wavefunction in Eq. ([5]) is in the 
instant form dynamics with equal- i in an ordinary frame, 
however the parton distributions are usually expressed 
in the infinite momentum frame or in the light-front dy- 
namics on the light-cone. Thus we need to transform the 
wavefunction in the instant form to that in the light-cone 
form. The connection between the two forms is through 
a Melosh-Wigner rotation (4(||43]. For a spin-^ particle, 
we use qr and gf to denote the instant and light-cone 
spinors, respectively, and the Melosh-Wigner rotation is 
known to be S3 



t 
q T 

i 
Qt 



p + + m 



P 



p + + m 



t 

I 
If 



m i/2 ( If 



t 
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If 



(6) 



where p R < L = p 1 ± ip 2 , p + = xAijj, ui — 

[(xM D + m) 2 + p 2 T ] ~ 1/2 with M 2 D = + 4±&, 

m and mo are mass parameters for the quark and di- 
quark, and M 1 / 2 denotes the Melosh-Wigner rotation 
matrix for the spin-1/2 particles. For the spin-0 scalar 
diquark, there is no such a transformation. For the spin- 
I vector diquarks, the transformation is represented by a 
3x3 matrix M l , whose explicit expression can also be 
found in Refs. [H, EH- In practice, we will not use the 
explicit form, for it does not appear in the final expres- 
sion for the spectator debris due to the unitary property 
M^M 1 = 1. 

Substituting Eq. © into §5§ and using Eq. ([3]), we 
obtain the TMD helicity distributions for the valence u 
and d quarks 



16tt 3 



cos 2 9tpgWs(x,PT) 



~7TH x \ siv ? e Vv w v(x,p T ), 
48-7T 3 3 

' "-' - ^3 \ sin2 9<PvW v (x,pt), 



gfl(x,p 2 T ) = - 



(7) 

(8) 
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where Wd is the Melosh-Wigner rotation factor for the 
helicity distribution with the form [t| |47| 



V-'H 



W D (x,p T ) 



(xM d + m) 2 - p 2 T 
(xM D + m) 2 + p\ 



2 • 



(9) 



Similarly the TMD unpolarized distributions are calcu- 
lated as 



^x(i sin 2 9<p 2 v + cos 2 0<p% ), (10) 



1 



1 



(11) 



The TMD helicity distributions thus are obtained by 
combining Eqs. ©, © QUI) and (fl~T|) . 



9iK x ,Pt) 



gtl{x,p 2 T ) = \f[ dv \x,p 2 T )W v {x, PT ). 



W s (x,p T ) 
(12) 
(13) 



Performing the integration over px, one arrives at the 
integrated helicity distributions 



9 r(x) = J d 2 p Tg r L (x,p 2 T ), 

gf{x) = J d 2 p T gH{x,p 2 T ), 



(14) 
(15) 



which have already been calculated in Ref. [9| in the light- 
con diquark model. 

Eqs. (IT2l and (|T3)) show that the helicity distribu- 
tions are connected to the unpolarized distributions in 
the light-cone diquark model. Here, we employ two ap- 
proaches to calculate g\h{x^p\\ corresponding two dif- 
ferent choice for f^(x,p^). In Approach 1 we adopt 
the Brodsky-Huang-Lepage prescription EH for the 
wavefunction in the momentum space 



if D = A D exp 




, (16) 



to calculate the unpolarized distributions flix^p 2 ^). In 
Approach 2 we factorize the ir-dependence and the px~ 
dependence of ff{x,Px) and choose the Gaussian form 
for the latter one in further 



f q 1 (x 7 p 2 T ) = f?(x) 



1 

/ 2 \q ex P 



-Pt 

(p 2 T ) q f 



(17) 



which has been used in many phenomenological appli- 
cations. For the integrated distribution (x) appearing 
in Eq. (|17p . we use the phenomenological input, i.e., the 
leading order set of the MSTW2008 parametrization [52j . 
In both approaches the pT-dependence for the polarized 
quarks is different from the Gaussian distribution which 
is adopted in Ref. 




FIG. 1: The 3-dimensional demonstration of xg^1{x,p^) (up- 
per panel) and — giL,{x,px) (lower panel) calculated by the 
light-cone diquark model in Approach 2. 



As a demonstration, in the upper panel and the lower 
panel of Fig. [IJ we show the 3-dimensional shapes of 
x 9il{ x >Pt) an( l ~ x 9il( x 'Pt) calculated from the light- 
cone diquark model using Approach 2, respectively. In 
this calculation the unpolarized valence quark PDFs 
fri x ) and ff v {x) are chosen at the scale Q = 1.2 GeV, 
and the Gaussian width for the unpolarized u and d 
quarks is adopted as (p 2 -)^ = 0.25 GeV f53j. Also we 
choose the mass parameters as m q = 0.35 GeV, Ms = 
0.6 GeV and My = 0.8 GeV, so that the sum of the quark 
and the diquark masses is larger than the proton mass in 
order to form stable bound state. 

Keeping the transverse momentum unintegrated, we 
can also explore the properties of the TMD distributions 
of quarks with positive helicity (denoted by q + (x,p 2 ^)) 
and those with negative helicity (denoted by q~(x,p^)), 
which can be obtained from the combination of ff(x, p\) 
and g\ L {x,p 2 T ), 



q f {x lP 2 T ) 
q i {x,p 2 T ) 



(f!(x,p 2 T ) + gf L (x,p 2 T )) , (18) 
(f!(x,p 2 T )-gl L (x,p 2 T )). (19) 



In Fig. [3J we plot the ratios (x , p 2 -,) j q^ (x , p\) for 
q = u v and d v as functions of pt at x = 0.3 (thick lines) 
and x = 0.15 (thin lines) in the light-cone diquark model 
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u T (x,p>-(x,p^) 



d , (x,p^)/d-(x,p^) 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
p T (GeV) p T (GeV) 

FIG. 2: The p T -profiles of ul(x,Pr) /ui(a;,pr) (left panel) 
and dl(x,p^)/di(x,p^) at x — 0.3 (thick lines) and x = 0.15 
(thin lines) in Approach 1 (solid line)and Approach 2 (dashed 
line) . 



using Approach 1 and Approach 2. From our model 
calculation we see that at small pr, the density of the 
valence u quark with positive helicity is several times 
larger than that with negative helicity at small px] but 
the difference between u\{x,p T ) and u^{x,p^) reduces 
with increasing px- For the d valence quark, the ratios 
of the positive helicity density and negative helicity den- 
sity are the same in both approaches. At smaller x and 
larger transverse momentum, the u quark has net nega- 
tive helicity, while the d quark tends to have net positive 
helicity. These features are originated from the quark 
transverse motion, which can be seen from the Melosh- 
Wigner rotation factor in Eq. (|9]) that is negative for large 
Pt and small x. 

In principle the TMD helicity distributions can be ac- 
cessed in SIDIS process (54| by employing the longitu- 
dinally polarized lepton beam off the longitudinally po- 
larized nucleon target, when the transverse momentum 
of the final fragmenting hadron is measured. The CLAS 
Collaboration at JLab has performed such measurement 
for the 7r + , 7T~ and 7r° productions [24]]. In the following 
we will use our model resulting TMD helicity distribu- 
tions to study the phenomenology of pion production in 
the double polarized SIDIS process. 

The differential cross section of the double longitudi- 
nally polarized SIDIS has the general form [5^, l56j 



da 



2ira 2 7 2 



dxdy dzdP? ± d(j) xyQ 2 V 2x/ ^ uu,t 
+ X e S ll B(e) F LL ), 



(20) 



where the higher-twist contribution has not been 
considered. The factors A(e) and B(e) are the depo- 
larization factors with the forms 



A(s) 



V 



2(1 -e) 



(21) 



B(e) = 



2(1 -e) 



Vl-e^y{l-±y), (22) 



The approximations hold in the limit 7 = 2Mx/Q — > 0. 
In the leading order of TMD factorization |57j the two 
structure functions Fjjut and Fll are expressed as 



Fi 



uu.t = x 



d 2 p T I d z K T S 2 (zp T - P h± + K T ) 



xf?(x,p T )Dl(z,K T ), 
F LL =xY / 4 Jd 2 p T J d 2 K T S 2 (zp T -P l 



h± 



X 



gl L (x,p 2 T )D\{z,K 2 T ). 



(23) 
K T ) 

(24) 



Here, Kt is the transverse momentum of the final pion 
with respect to the momentum of the fragmenting quark, 
and Ph± is the transverse momentum of the pion with re- 
spect to the photon momentum in the target rest frame. 
For the TMD FFs needed in the calculation, and we as- 
sume that their i^-dependence has a Gaussian form 



D\ (z,K 2 )=Dl(z) 



1 



tt(K, 



^eM~K 2 /(K 2 )), (25) 



and we choose the Gaussian width (Kj.) — 0.20 GeV [53| 
for the u and d quarks. For the collinear FFs we adopt 
the leading-order set of the Kretzer parametrization [58| . 
The DSA in SIDIS can be measured through 



A L l(x, z, Pt) 



da n - d<T n - da^ + da u 
B(e)F LL 



da^ 
1 



xyQ 2 



dan 

1 + ; 



xyQ 



1 fuf 

W V 1 + 2x 



A{e)F l 



(26) 



UU,T 



where f/i denotes the positive/negative helicity of the 
lepton and the nucleon, respectively. On the other hand, 
one may define the DSA that is scaled by the depolariza- 
tion factors 



A\ 



A(s) 



b{£)All ^(i+fi)W 



xyQ 



1 



2.T 



LL 



(27) 



which is in consistency with the measurement by the 
CLAS Collaboration 

In Fig. [3l we plot the P/j_L-dependent DSAs defined in 
Eq. (|2T|) for the 7r + , ir~ and ir° productions at the CLAS 
kinematics 

£ e = 5.7GcV, 0.12 < x < 0.48, y < 0.85, 

0.4 < z < 0.7, 0.9 GeV 2 < Q 2 < 5.4 GeV 2 , (28) 

W 2 > 4 GeV 2 , M x (eTr°) > 1.4 GeV. 



and compare them with the CLAS data [24{ . The dashed 
and the solid curves correspond to the predictions using 
the TMD unpolarized and helicity distributions from Ap- 
proach 1 and 2 of the light-cone diquark model, respec- 
tively. It is found that the TMD distributions in Ap- 
proach 2 agree with the CLAS data generally, that is, we 
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1 i I i I i I i I r , | i | i | , | J!-, L 



~ i — 1 — i — 1 — r 




-0.1 

0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

P hi (GeV) 

FIG. 3: Double spin asymmetries A\ for the tv + , tt~ and ty° 
productions as functions of the pion transverse momentum 
Phx ■ The dashed and the solid curves correspond to the pre- 
dictions using the TMD unpolarized and helicity distributions 
from Approach 1 and 2 of the light-cone diquark model, re- 
spectively. Data are from [24[. 



Ref. 37]. Boer, Gamberg, Musch and Prokudin devel- 
oped a framework to avoid the convolution of the trans- 
verse momenta through the so-called Bessel-weightcd 
procedure. In this framework, instead the TMD dis- 
tributions, the Fourier transformation of them (denoted 
by f q (x,b^)) can be accessed directly by measuring the 
Bessel-weighted asymmetries, where bx is the transverse 
position of partons conjugate to pr- Thus f q (x, b 2 -,) are 
also referred to as the parton distributions in the trans- 
verse position space. Once f q (x,b^) are measured, the 
TMD distributions can be obtained by the inverse Fourier 
transformation on f q (x,b^)- The Fourier-Bessel trans- 
forms is also useful in studying the factorization as well 
as the scale dependence of TMD cross section (59l - [65j . 

In this section, we study the helicity distribution in the 
transverse position space and the related phenomenology 
in SIDIS processes using the Bessel-weighted framework 
developed in Ref. [13] ■ Using the identity 





(p 2 t)7 


<Pr)f 


(p 2 t )7 




Approach 1 


0.25 GeV 2 


0.25 GeV 2 


0.14 GeV 2 


0.13 GeV 2 


Approach 2 


0.10 GeV 2 


0.10 GeV 2 


0.08 GeV 2 


0.08 GeV 2 



TABLE I: The mean square momenta of the unpolarized and 
the polarized quarks in Approach 1 and 2. 



predict that the asymmetries for all three pion channels 
have the tendency to fall off at large Phx- However, the 
model has the difficulty to explain the positive slope for 



7r~ production at moderate Pt region. In Ref. [24[ the 
Py-dependence of the double-spin asymmetry was inter- 
preted by difference widths of the quark TMD distribu- 
tions with different flavor and polarizations [III ] resulting 
from different orbital motion of quarks polarized in the 
direction of the proton spin and opposite to it [l2| . In 
Table. Inl we list the widths for the polarized quarks that 
are defined as 



1x17 dx [ d2pT p T 91l <>> Pt) 
XTrr dx I d2 PT9 q iL 0: Pt) 



(29) 



in the light-cone diquark, together with the widths for 
the unpolarized quarks. The integrating limits for x are 
chosen as x m i n = 0.1 and x max — 0.9, corresponding 
to the valence region. In both cases, the widths for the 
polarized u and d valence quarks are almost the same. 
The ratio (p T ) g /(p T ) / is around 0.52 — 0.56 in Approach 
2, while that in Approach 1 is 0.8. 



III. HELICITY DISTRIBUTION IN FOURIER 
SPACE AND BESSEL-WEIGHTED DSA 



Because the structure functions given in Eqs. ([23|) and 
involves the convolution of the transverse momenta 
of partons, it is not easy to obtain directly the informa- 
tion of TMD distributions when P^± is measured. In 



S 2 (z Pl 



h± 



K T ) 



(2tt 



"£_ e %bT-{zpT— PhA 



-K T ) 

(30) 



the integrations over the transverse momenta in Eqs. ([23 
and (EU) are deconvoluted [H 



UU,T 



\b T \J (\b T \\P, 



h± 



xf?(x,z 2 b T )Dl(z,b 2 T ), 
d\b T ' 



2n 



b T \J (\b T \\P hl _\) 



xg q L (x,z'b T )Dl(z,b T ), 



(31) 



(32) 



where the functions with a tilde are the Fourier trans- 
forms of the TMD distributions and FFs 

~ft{x : z 2 b 2 T ) = j ' d 2 p T e izb ^ fl{x,p T l (33) 

~gl L (x,z 2 b 2 T ) = J d 2 p T e lzb ^gl L (x,p 2 T ), (34) 

Dl(z,b 2 T ) = J d 2 p T e lb ^D\{z lP 2 T ). (35) 

The Bessel-weighted asymmetry corresponding to the 
DSA defined in Eq. (|2T|) can be obtained by the weighting 
function J (|#T||-P/ii|), w ^ &T the transverse coordi- 
nate conjugate to Ph± [37| 

A M\B T \\P h ^) {BT) 

^(l + £)fd\P h± \\P hX \J (\B T \\P h± \)F LL 



xyQ 
1 

tiQ 



l+£)f d \P> 



h± 



\M\BT\\Ph±\)F UUt T 



^(l + £) E q e 2 a ~9 q 1L ^^T)Dl(z,B 2 T ) 



v&{ 1 + £) E q e 2 fKx,^ 2 T )Dl( Z) B 2 T ) 



(36) 
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where convolutions of TMD PDFs and FFs become sim- 
ple products. Another advantage of the Bessel-weightcd 
asymmetries is that the universal soft factors appearing 
in the numerator and the denominator of Eq. (|36[) can- 
cel [I?} ■ In principle the hard parts Htl and Huu,t 
for the polarized and unpolarized subprocess should 
also present in the numerator and the denominator of 
Eq. (|36p . At higher order they may receive different cor- 
rections. However, the calculation of the hard parts at 
one-loop order shows that Hll = Hjju.t- Thus it is 
a quite good approximation to use Eq. (|36[) to perform 
numerical calculation. Furthermore, in the case of Ph±- 
dependent asymmetry the leading-order expression (|27|) 
can describe the data, as we have shown in the previous 
section. 

We use Approach 2 in the light-cone diquark model to 
study the Bessel-weighted DSA. According to Eqs. (|39|) 
and (|25[) . in this approach the unpolarized distribution 
and FFs in the Fourier space have the analytic forms 



d 2 PT j bTPT 



fl(x) exp 



f_M 



exp 



-Pt 



(37) 



D\(z,b T ) = / d z K T e 



D\{z) exp 



-b 2 T (K* T ) 




FIG. 4: The 3-dimensional demonstration of xg"i,(a 



(up- 



per panel) and — g\t{x, b T ) (lower panel) calculated by the 
(38) light-cone diquark model in Approach 2. 



Thus the Gaussian widths for f1(x,b^) and D\(z,b^) 
are ^/{p T ) q j and 4/(if|,), respectively. In our model the 
helicity distributions in the Fourier space can not be ob- 
tained analytically. We calculate them numerically in- 
stead, and show the 3-dimensional shapes of (a;, b\) 
and —xgf" L {x,b\) in Fig. |U We see that both of them 
are sizable in the region br < 1 fm, and approach to zero 
at bx = 2 fm. 

The Bessel-weighted DS As as functions of Bt for tt + , 
ir~ and 7T° productions at the CLAS kinematics (|2"8)> 
are plotted in Fig. [SJ The solid curves are the results 
from Approach 2 in the light-cone diquark model. As a 
comparison, we also present the prediction by using the 
TMD helicity distributions adopted in Ref. 
gi(x,p^) has a Gaussian dependence for px 



111, where 



9il(x,Pt) = 9%( x ) 



1 



t (Pt 



rg-exp 

19 



-Pt 



(39) 



and g\{x) is taking from the GRSV2000 leading-order 
parametrization Q. The results are show by the dashed 
and dotted curves for (p|)| = 0.17 GeV 2 and 0.10 GeV 2 , 
respectively. Although the TMD helicity distributions 
from the light-cone diquark model (Approach 2) and from 
the Gaussian ansatz can both describe the CLAS data, 
their predictions for the Bessel- Weighted DSAs are ap- 
parently different, especially for the 7r + production. Thus 
we expect that the measurement of the Bessel-weighted 
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FIG. 5: The Bessel-weighted DSAs A( o(]BTllPh±l) for tt + , 
7r~ and 7T° productions as functions of Bt at CLAS. The 
solid lines are from the Approach 2 of the light-con diquark 
model, while the dashed line and the dotted lines are from 
the Gaussian ansatz for the TMD helicity distributions with 
(Pt)1 = 0.17 GeV and 0.10 GeV 2 , respectively. 



DSAs can shed light on the helicity distribution in the 
Fourier space. 

In calculating the curves in Fig. El we have used the last 
line of Eq. (|36[) where the whole range of Ph± has been in- 
tegrated over. As a check, we also calculate ^°^ BT ^ Ph± ^ 
using the integration limit GeV < Ph± < 1.12 GeV ac- 
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FIG. 6: The sizes of the numerator (solid line) and the de- 
nominator (dashed line) of Eq. (|36[) as functions of Bt for 7r + 
production. 



cording to the kinematical cuts at CLAS. We find that 
the difference of these two results is about one percent, 
which agrees with the conclusion in Ref. [37j that the 
contribution of the large Ph± tail is suppressed in the 
Bessel-weighted integrals. This is important such that 
Eq. (j3"6")l can be used in phenomenological analysis to 
extract giL{x,q T ) directly without worrying about the 
large Phi. contribution. 

In Fig. [5J we plot the sizes of the numerator and the 
denominator of Eq. (l36l) as functions of Bt for 7r + pro- 
duction, calculated from the light-cone diquark model. 
The result shows that the measured differential cross sec- 
tion in the Fourier space will decrease almost exponen- 
tially with increasing Bt- Thus the region where Bt is 
not so large can be explored in SIDIS. In the case the 
Bt < 1.5 fm region is measured, it will provide the infor- 
mation of giL(x, bj.) in the region &x = zBt < 1 fm at 
CLAS, where giL(%, b T ) should be sizable. 

An interesting phenomenon is that all the Bessel- 
weighted asymmetries for the three pions increase with 
increasing Bt- This can be explained by the fact that 
the mean square of &r for the polarized distribution in 
the Fourier space 



2)q = SxZ" dx f d2bT b T 9il ( x > b r) 

9 s::::dxs<pb T ~gi L {x,b%) ■ 



(40) 



is larger than the width (b T ) 9 ^ of the unpolarized distribu- 
tion, which is opposite to the case of TMD distributions 
because 6t is conjugate to pr- Explicitly, in Approach 
2 of the light-cone diquark model we obtain (b T )g V = 
1.45 fm 2 , (b T )™ = 1.40 fm 2 , and (b T )j = 0.62 fm 2 . 

DSAs for meson production in SIDIS as functions of 
x have also been measured at HERMES [38[ by scat- 
tering a 27.6 GeV positron beam off longitudinally po- 
larized hydrogen and deuterium targets, as well as at 
COMPASS [a3 by collision of a 160 GeV muon beam 
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FIG. 7: Same as Fig. but for the HERMES kinematics. 
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FIG. 8: Same as Fig. [5] but for the COMPASS kinemati 



on a longitudinally polarized NH 3 target. In Fig. [7] and 
Fig. [51 we present the prediction for the Bessel-weighted 
DSAs of 7r + , 7r~ and tt° productions off a proton target 
in SIDIS at HERMES and COMPASS, using Approach 
2 of the light-cone diquark model as well as the TMD 
helicity distributions from Ref. [ll| for comparison. The 
kinematical cuts applied in the calculations are 

0.023 < x < 0.6, y < 0.85, 

0.2 < z < 0.8, Q 2 > 1 GeV 2 , W 2 > 10 GeV 2 (41) 

for HERMES 0, and 

0.004 < x < 0.7, 0.1 <y< 0.85, 

0.02 < z < 1, Q 2 > 1 GeV 2 , W 2 > 25 GeV 2 (42) 

for COMPASS [Hj|. Again we observe that the Bessel- 
weighted asymmetries for the three pions increase with 
increasing B T at HERMES and COMPASS. The differ- 
ence is that the asymmetries at HERMES and COM- 
PASS are smaller than that at CLAS, this is because the 
HERMES and COMPASS can probe the smaller x re- 
gion where the helicity distributions are small compared 
to the unpolarized TMDs. 
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IV. CONCLUSION 

In summary, we study the helicity distributions of va- 
lence quarks in the transverse momentum space and the 
transverse coordinate space. The TMD helicity distri- 
butions of valence u and d quarks are calculated in the 
light-cone diquark model by adopting two different ap- 
proaches. By comparing the pT-dependence of the he- 
licity distributions at different x, we find that at small 
x and large pr, the valence u quark has net negative 
helicity while the d quark tends to have net positive he- 
licity. We thus use our model results to analyze the Ph±- 
dependent DSAs of 7r+, tt~ and 7r° productions in SIDIS, 
and compare them with the data measured by the CLAS 
Collaboration. We find that the asymmetries calculated 
in Approach 2 of the light-cone diquark model agree with 
the CLAS data. We then calculate the valence helicity 
distributions in the Fourier space using Approach 2 of 
the light-cone diquark model, and use them to predict 
the Z?T-shape of Bessel- weighted DSAs for ir + , 7r~ and 
7T° productions in SIDIS at CLAS, COMPASS and HER- 
MES for the first time. Because the mean square trans- 



verse position of the polarized quarks is larger than that 
of the unpolarized quarks, the Bessel-weighted asym- 
metries for all three pions increase with increasing Bt- 
The contribution from large Ph± region to the Bessel- 
weighted DSAs are calculated numerically and is found 
to be negligible. Our study suggests that the Bessel- 
weighted asymmetry ^ ^ BT ^ F ' h± ^ mav be measured in 
the double polarized SIDIS and provides a straightfor- 
ward access to the helicity distributions in the coordinate 
space. 
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